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Abstract 

We resolve the B — ► <pK* polarization puzzle by postulating a smaller B — ► K* form 
factor Aq ~ 0.3 and by adding penguin annihilation and nonfactorizable contributions from 
the perturbative QCD approach. If this explanation is valid, the penguin-dominated modes 
governed by the B — ► K* form factors, such as B + — > an d 5° — > ET* ^* , should 

exhibit similar polarization fractions. Our resolution is compared with others in the literature, 
and experimental discrimination is proposed. 
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To understand the polarization fractions of the B — > <j)K* decays has been a challenge. Motivated 
by this subject, we have investigated most of the B — > VV modes, and observed that they are 
classified into four categories [1]. First, the B° — > (D* s + , D* + , p + )D*~ modes can be understood by 
kinematics in the heavy-quark limit, whose longitudinal polarization fractions Rl ~ 0.5, 0.5, and 
0.9 [2, 3], respectively, follow the mass hierarchy among the D*, D* and p mesons emitted from 
the weak vertex. Second, the B — > (p,u)p modes are understood by kinematics in the large-energy 
limit, leading to R L ~ 1 [4, 5, 6]. Applying the same estimation, we have predicted R L ~ 0.7 for 
the B + — > (D* + , D* + )p° decays, which can be compared with future data. For penguin-dominated 
modes, such as those listed in Table 1, the polarization fractions deviate from the naive counting 
rules based on kinematics [7]: the annihilation contribution from the (S — P)(S + P) operators and 
the nonfactorizable contribution decrease Rl to about 0.75 for the pure-penguin B + — > p + K*° decay. 
Adding a tree contribution, R L of B + — > p°K* + can go up to about 0.9 [1]. The fourth category, 
consisting of the puzzling B — > <J)K* decays, is also pure-penguin, but its R L ~ 0.5 shown in Table 1 
is much lower than 0.75. 





Mode 


Pol. Fraction 


Belle 


Babar 


B + 


-> <pK* + 


Rl 

R± 


0.49 ±0.13 ±0.05 [8] 
0.1218& ±0.03 [8] 


0.46 ±0.12 ±0.03 [9] 


B° 


-> (j)K*° 


Rl 


0.52 ±0.07 ±0.05 [8] 


0.52 ±0.05 ±0.02 [10] 






R± 


0.30 ±0.07 ±0.03 [8] 


0.22 ±0.05 ±0.02 [10] 




Mode 


Pol. Fraction 


Belle 


Babar 


B+ 
B + 


-> p°K* + 
-> p+K*° 


Rl 
Rl 


0.50 ±0.19lg;g? [11] 


0.96t£S£ ± 0.04 [9] 
0.79 ±0.08 ±0.04 ±0.02 [12] 



Table 1: Polarization fractions in the penguin-dominated B — > VV decays. 

It seems that the 5 — > polarizations are the only anomaly so far, and many attempts to 
resolve it have been proposed, which include new physics [13, 14], the annihilation contribution [15] in 
the framework of QCD-improved factorization (QCDF) [16], the charming penguin in soft-collinear 
effective theory (SCET) [17], the rescattering effect [18, 19, 20], and the b — > sg transition (the 
magnetic penguin) [21]. We have carefully analyzed these proposals [1]: the annihilation amplitude 
has to be parameterized in QCDF, and varying free parameters to fit the data can not be conclusive 
[22]. The charming penguin strategy, demanding many free parameters, does not help understand 
dynamics. Moreover, it has been argued that the charming penguin, without infrared divergences 
from diagrammatic calculations, should be factorizable in the current leading-power SCET formalism 
[1]. A similar criticism has been raised recently in [23]. The rescattering effect is based on a model- 
dependent analysis [24, 25], and constrained by the B — > pK* data. The prediction R» ^> R± for 
B — > <f)K* [20], R\\ and R± being the parallel and perpendicular polarization fractions, respectively, 
also contradicts the observed relation R\\ ps R^ in Table 1. The exotic magnetic penguin is suppressed 
by the G-parity, and not sufficient to reduce Rl down to 0.5 [1]. However, we are not claiming a 
signal of new physics, since the complicated QCD dynamics in B — > VV decays has not yet been 
fully explored. 
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In this letter we shall investigate whether QCD effects can resolve the B — * <pK* polarization 
puzzle without resorting to exotic mechanism or new physics. These decays have been studied in 
the perturbative QCD (PQCD) approach [26, 27, 28], and the results of the branching ratios, the 
magnitudes of the helicity amplitudes A L , A\\, and A±, and their relative strong phases 4>\\ an d <f>± 
are summarized in Table 2 [7]. The normalization of these amplitudes have been chosen, such that 
they satisfy 

\A L \ 2 + \A l{ \ 2 + \A ± \ 2 = l, (1) 

with \A L \ 2 = R L , \A\\\ 2 = and \A±\ 2 = R±. The first rows (I), coming only from the factorizable 
emission topology, correspond to the results under the factorization assumption (FA) [29]. It is 
obvious that the polarization fractions Rl ~ 0.92 and R\\ R± ss 0.04 follow the naive counting 
rules, 

R L ~1- 0{m\lm\) , R\\ ~ R ± ~ 0(mj/m|) , (2) 
m b (m^) being the B (0) meson mass. 



Mode 


Br (10~ y ) 


\Al\ 2 


l^lll 2 


|4l| 2 


(f>\\(rad.) 


(f)±(rad.) 


<f)K* (l) 


14.48 


0.923 


0.040 


0.035 


7T 


7T 


(II) 


13.25 


0.860 


0.072 


0.063 


3.30 


3.33 


(HI) 


16.80 


0.833 


0.089 


0.078 


2.37 


2.34 


(IV) 


14.86 


0.750 


0.135 


0.115 


2.55 


2.54 


<f>K* + (l) 


15.45 


0.923 


0.040 


0.035 


7T 


7T 


(II) 


14.17 


0.860 


0.072 


0.063 


3.30 


3.33 


(HI) 


17.98 


0.830 


0.094 


0.075 


2.37 


2.34 


(IV) 


15.96 


0.748 


0.133 


0.111 


2.55 


2.54 


<j)K* 


10.2t£? 


u - oy -0.02 


n 22^ U - Ui 

u - zz +0.01 


u - iy +0.01 


o 09+O.II 


9 o-,+u.r2 



Table 2: (I) Without the nonfactorizable and annihilation contributions, (II) add only the nonfactoriz- 
able contribution, (III) add only the annihilation contribution, and (IV) add both the nonfactorizable 
and annihilation contributions. The last row is for A = 0.28. 

The next-to-leading-power annihilation amplitudes, mainly from the (S — P)(S + P) operators, 
and the nonfactorizable amplitudes bring the first rows into the fourth ones (IV) with the fractions 
Rl ~ 0.75. We observe from the second and third rows, (II) and (III), that these subleading 
corrections work toward the direction indicated by the data. It is easy to understand the sizable 
deviation from Eq. (2) caused by these subleading corrections, which are of Oim^/mB) for all the 
three final helicity states [7] . If they are of the same order of magnitude as and constructive to the 
transverse polarization amplitudes, an enhancing factor will be gained, which may be large enough 
to modify the counting rules numerically (note that m^/mB is only about 1/5). However, the total 
effect, as shown in Table 2, is not sufficient to lower Rl of the B — > <pK* decays down to around 0.5. 
The branching ratios in (I) and in (IV) are roughly equal, indicating that the subleading corrections 
decrease the longitudinal components and increase the transverse ones by roughly equal amount. 
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Two nice features exhibited in Table 2 are that PQCD has predicted R\\ ~ R±, contrary to those 
from the rescattering effect [20], and that the relative strong phases among the helicity amplitudes 
are consistent with the B° — > (f)K*° data: 

0„ = 2.21 ± 0.22 ± 0.05 (rad.) , <f> ± = 2.42 ± 0.21 ± 0.06 (rad) [8] , 

011 = 2.34±g;|g ± 0.05 (rad.) , <f> ± = 2.47 ± 0.25 ± 0.05 (rad) [10] . (3) 

The former implies that the rescattering effect may not be essential in B meson decays into two light 
mesons [30]. The consistency of the predicted 0m and <f)± with the data, once again, supports that 
the evaluation of strong phases in PQCD is reliable. Other examples include the predictions for the 
direct CP asymmetries in the B — > K + ir~ , it + it~ modes [27, 28], and the results of the B — > D^n(p) 
branching ratios, which crucially depend on the strong phases of the color-suppressed amplitudes. 

As emphasized above, the B — > <pK* polarizations are very unique, and it is difficult to find new 
mechanism, which affects only these modes but not others. Hence, we do not intend to propose any 
new mechanism or new physics to resolve the puzzle. To explain our idea, we quote the explicit 
expressions of the three helicity amplitudes in terms of the B — > K* transition form factors in FA 

[1], 

A L oc 2r 2 e* 2 (L) ■ e* 3 (L)A , (4) 
A\\ oc -V2(l + r 2 )A 1 , (5) 

A ± oc -^^2{(v 2 .v^-l]V, (6) 

with the K* (<fi) meson velocity v 2 (v 3 ) and polarization vector e 2 (e^), r 2 = mx*/rriB and r% = 
m^/mB- The form factors Aq, A±, and V in the standard definitions obey the symmetry relations in 
the large-energy limit [31, 32], 

m^ v , = = = ^ 

m B + m K * 2E 2E 

m K * m B + m K * m B - m K * 

~^ Ao ~ 2E Al ^b~ A2 ' (8) 

where T\ and T 2 are the form factors involved in the B — > K*^ decays, and E is the K* meson 
energy. 

The results in Table 2 correspond to the form factors A = 0.40, A 1 = 0.26 and V = 0.35. First, 
the B — ■> K*^( branching ratios have constrained the form factors 7\ ~ T 2 pa 0.3 [33, 34], which are 
also in agreement with the lattice result [35]. Compared to the symmetry relation in Eq. (7), it is 
obvious that PQCD has given reasonable values of A 1 and V. Second, there has not yet been any 
measurement, except B — >■ <f)K*, which constrains Aq. The other penguin-dominated B — > p(u)K* 
decays are mainly governed by the B — > p(u) form factors. Third, the PQCD predictions for the 
B — > (j)K* branching ratios in Table 2 are larger than the data [36], 

B(B° ^(f)K*°) = (9.5 ±0.9) x 10~ 6 , B(B + -> <f)K* + ) = (9.7 ± 1.5) x 10" 6 . (9) 

Note that the same value of Aq pa 0.40 leads to the branching ratio about 10 x 10 -6 for the longitudinal 
component in PQCD, but about 5 x 10 -6 in QCDF [15, 37], because of the dynamical penguin 
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enhancement in the former [27]. The above three observations hint that the PQCD results for 
the transverse components of the B — > 4>K* decays should have been reasonable, and that the 
longitudinal components may have been overestimated. We are then led to conjecture that a smaller 
A will resolve the puzzle, giving both lower Rl and lower branching ratios. 

In PQCD, a B — > K* form factor is written as the convolution of a hard kernel with the B meson 
wave function and with a set of K* meson distribution amplitudes. Note that the form factors A , 
A 1 and V involve different sets of K* meson distribution amplitudes: the twist- 2 <p K * , and the two- 
parton twist-3 <ft K » and <j) s K * for A , and the twist- 2 and the two-parton twist-3 (f> v K * and for 
A\ and V (the notations are referred to [7]). Our investigation indicates that the latter set of model 
distribution amplitudes derived from QCD sum rules [38] has been acceptable, but the former set 
has not. Recently, the reanalysis of (fix*, parameterized as 

<f>K*(x) = ^^x(l-x)[l + 3af(l-2x)] , (10) 

showed that the Gegenbauer coefficient has a revised value af* = 0.10 ± 0.07 [39], different from 
af* = 0.19 ± 0.05 in [38]. That is, considering the theoretical uncertainty, 4>k*{x) could be quite 
close to the asymptotic model corresponding to af* = 0. 

To test our idea, we choose the asymptotic models for the K* meson distribution amplitudes 
relevant to the evaluation of A : 

<l>K'{x) = ^^X(l-X), (11) 

= tMw^-Z*) 2 , (12) 



2y/W c 

Ik* 



i>' K .(x) = -%=3(l-2x), (13) 



which lead to A = 0.28, about 70% of the original value. The main reduction is caused by the change 
of 4> S K * in Eq. (13). The model-dependent evaluations of A vary in a wide range from 0.31 to 0.47, 
and A rs 0.3 has been supported by the recent covariant light-front QCD (LFQCD) calculation [40]. 
This smaller value does not contradict to any existing data as emphasized above. We suggest to also 
reanalyze <$ K * and <f> s K * in QCD sum rules, so that it is possible to examine whether a consistency 
between the PQCD and LFQCD results of A can be achieved. 

The models for the distribution amplitudes 4>k*i $"k* an d 4>°k*i relevant to the evaluation of the 
form factors A 1 and V, and those for the meson distribution amplitudes and for the B meson wave 
function, remain the same as in [7]. We then compute all amplitudes, including the penguin emission, 
penguin annihilation and nonfactorizable ones, for the longitudinal and transverse polarizations using 
the kr factorization formulas in [7]. The numerical outcomes are listed as the last row in Table 2. 
Simply adopting the asymptotic models in Eqs. (11)-(13), the modified branching ratio 10.2 x 10~ 6 , 
the polarization fractions R L = 0.59 and R\\ xs R±, and the relative strong phases 4>\\ ~ 4>± ~ 2.3, are 
all consistent with the B° — > (pK*° data in Table 1, and in Eqs. (3) and (9). Therefore, we claim that 
the measured B — > <pK* polarizations might imply nothing but a smaller form factor Aq, and that 
their explanation does not require any exotic mechanism or new physics. The penguin annihilation 
and nonfactorizable contributions play an important role here. In FA without these contributions, 
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Aq has to be as small as 0.15 in order to reach Rl ~ 0.6, for which the B — > <pK* branching ratios 
will fall far below the data. 

The central values in the last row of Table 2 correspond to the shape parameter uj b = 0.40 GeV, 
appearing in the B meson wave function [27], 



and the variable b conjugate to the parton transverse momentum in the B meson. The errors in 
the superscripts (subscripts) come from oj b = 0.36 {uj b = 0.44) GeV. The range uj b = 0.40 =F 0.04 
GeV, determined by a fit to the B — > n form factor from light-cone sum rules [41, 42], leads to the 
B -> K* form factors A = 0.28±g;g|, A 1 = 0.26^, and V = 0.35±g;g|. The above theoretical errors 
simply mean those arising from the unknown B meson wave function. It is clear that the polarization 
fractions are insensitive to this source of uncertainties. There are certainty other sources of theoretical 
uncertainties, whose detailed investigation is not the focus of this work. 

If our explanation is valid, the B — > <pK* decays can be classified into the same category as of 
B — > pK*, for which the sub leading penguin annihilation and nonfactorizable contributions render 
the polarization fractions deviate from the naive counting rules based on kinematics. The only 
difference is that the form factor ratio Aq/ A x for the B — > K* transition is smaller than that for the 
B — > p transition. For example, LFQCD gave A /Ai = 1.19 for the former, and A /Ai = 1.27 for the 
latter [40]. PQCD, using the set of p meson distribution amplitudes from [38], gave Aq/A\ = 1.63 
[41]. The PQCD results for the B — > np decays, to which the B — > p form factors are relevant, 
have been in good agreement with the data [43]. It implies that the values of the B — > p form 
factors yielded in PQCD are reasonable. Therefore, the longitudinal polarization fraction R L of the 
B + — »■ p + K*° decay should be larger than that of B — > <pK* . An explicit study of the B — > pK* 
modes in PQCD is in progress. Furthermore, if our explanation is correct, the modes governed by 
the B -> K* form factors, such as B+ -> K* + K*° and B° -> K*°K*°, must have low R L « 0.6. Note 
that the B — > <f>K* modes occur through the b — > s penguin, while the B —> K*K* modes occur 
through the b — > d penguin. From the viewpoint of PQCD, the B — > uoK* decays, without exotic 
mechanism, show R L close to those of B — > p JT*. All the above predictions can be confronted with 
the future polarization measurement of the B — > K*K*, uK* decays. 

At last, we compare the polarization fractions of the B — > pK*, ouK* decays derived from the 
different approaches with the B — > <f)K* data in Table 3. The prediction from the rescattering effect 
[18, 19, 20] is easily understood: the B pK*, <pK* decays involve the same £)(*)£)(*) intermediate 
states, and their polarization fractions are certainly almost equal. The relation R\\ (4>K*) 3> R±((f)K*) 
is attributed to the vanishing in the B — > D*D* channels. The b — > sg transition contributes 
to the B — > ujK*, <f)K* modes, such that their polarization fractions are close to each other [21]. 
Without the rescattering effect, the parallel and perpendicular polarization fractions are expected to 
be similar. This is also the case in PQCD [7] and in QCDF [15]. The future data can provide an 
unambiguous discrimination among these proposals. We did not list QCDF in Table 3, because both 




(14) 



where the normalization constant N B is related to the decay constant f B through 




(15) 
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the relations R L (p + K*°) > R L (4>K*°) and R L (p + K*°) R l (4>K*°) are allowed within the involved 
large theoretical uncertainty [44]. Hence, it is difficult to discriminate QCDF experimentally from 
the others. 



PQCD 


R L (pK*) 


> R L (<PK*) 


R\ 


\{4>K*) 








Rl{ooK*) 


» R L (p°K*) 










rescattering 


Rl(pK*) 




R\\ 


{<t>K*) 


> i?_ 


L(0^*) 


magnetic dipole 


Rl(ujK*) 




R\ 


\{4>K*) 




_(<^*) 



Table 3: Comparison of the B — > plT", cjif* polarization fractions from the different proposals with 
the B -> data. 

In this letter we have proposed a possible resolution to the B — > polarization puzzle within 
the Standard Model, which was found to be nothing but the consequence of a smaller B — > K* 
form factor A . Postulating A m 0.3, which does not contradict to any existing measurement, we 
are able to explain the data in the PQCD approach. This form factor value leads to Rl = 0.84 in 
FA (smaller than 0.92 corresponding to A = 0.40 in Table 2). The penguin annihilation from the 
(S — P)(S + P) operators and the nonfactorizable contribution, which can be estimated reliably in 
PQCD, then further bring R L down to 0.59. The remaining concern is whether the asymptotic models 
in Eqs. (11)-(13) are allowed within theoretical uncertainty. We have suggested to reanalyze these 
K* meson distribution amplitudes appearing in the factorization formula for Aq in the framework of 
QCD sum rules. Because of the unknown Aq, which is a source of QCD uncertainty, it is still too early 
to claim any exotic mechanism or new physics in the B — * <pK* polarization data. We have predicted 
that the penguin-dominated decays governed by the B — > K* form factors, such as B + — > K* + K*° 
and B° — > K*°K*°, should exhibit similar Rl ~ 0.6. The comparison among the different proposals 
for resolving the puzzle has been summarized in Table 3, which can be discriminated experimentally 
in the near future. 

I thank C.H. Chen for the numerical analysis, and H.Y. Cheng, C.K. Chua, A. Kagan, and S. 
Mishima for useful discussions. This work was supported by the National Science Council of R.O.C. 
under Grant No. NSC-93-2112-M-001-014 and by the Taipei Branch of the National Center for 
Theoretical Sciences of R.O.C. 
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